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ABSTRACT 

We present in this letter a high resolution (22" FWHM) extended map at 2.1mm of the 
Sunyaev-Zel'dovich effect toward the most luminous X-ray cluster, RXJ1347-1145. These obser- 
vations have been performed with the D1AB0L0 photometer working at the focus of the 30m 
IRAM radiotelescope. We have derived a projected gas mass of (l.l±0.1)xl0 14 /i 50 5/2 M Q within 
an angular radius of 9 = 74" (ie: projected radius of 0.6Mpc, Hq = 50km/s/Mpc, Sl m = 0.3, 
Ha = 0.7). This result matches very well the expected gas mass from the cluster models of X-ray 
data. With an unprecedented sensitivity level our measurement does not show significant depar- 
ture from a spherical distribution. The data analysis also allows us to characterize the 2.1mm 
flux of a well known radio source lying in the center of the cluster: _FRs(2.1mm) = 5.7 ± 1.6mJy. 

Subject headings: cosmology: cosmic microwave background — cosmology: observations — galaxies: 
clusters: individual (RXJ1347-1145) — intergalactic medium 

1. Introduction Sadat, Blanchard and Oukbir 1998; Bahcall and 

Fan 1998). This strong coupling makes the study 

The statistical properties of galaxy clusters of massive and distant dusters a very useful tool 

(shape, structures, size, temperature, mass) de- f or cosm ology. The intergalactic gas component 

pend strongly on the geometry of the Universe. can be observed through its Bremsstrahlung emis- 

Their study provides some robust constraints on sion at x _ ray wave l en gths. It can also be de- 

cosmological models, structure formations and tected from su bmillimeter to radio wavelengths 
evolution (Oukbir, Bartlett and Blanchard 1997; 
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via the Sunyaev-Zel'dovich effect (Sunyaev and 
Zel'dovich 1972). Whereas the X-ray emission 
depends on the square of the gas density, the 
Sunyaev-Zel'dovich (SZ hereafter) effect is linearly 
dependent on this quantity. For this reason, the 
SZ effect is proportional to the column density 
and thus to the line of sight integrated gas mass. 
Moreover the SZ signal does not suffer from the 
brightness decrease of the radiations (in partic- 
ular the X-ray one) due to the expansion of the 
Universe (see a review by Birkinshaw (1999)). 

The RXJ1347-1145 cluster is known as the most 
luminous X-ray cluster to date. It has been de- 
tected in the ROSAT All Sky Survey and fur- 
thermore studied with the ROSAT-HRI and the 
ASCA-GIS2 instruments (Schindler et al. 1995, 
1997). With an intrinsic bolometric X-ray lumi- 
nosity of L bo i = 21 x 10 45 h§ ergs s _1 , it shows a 
very peaked X-ray emission profile with an angu- 
lar core radius of 9 C = 8.4 ± 1.8". It also presents 
a very strong cooling flow in its central region, 
with a corresponding accreting rate of M coo j = 
3OOOM /yr. It is a distant cluster with z = 
0.45, so that the angular distance is 1670 /ig Mpc 
(490pc/arcmin, with H = 50km/s/Mpc, TO = 
0.3, f^A = 0.7; 410pc/arcmin in a standard CDM 
model). The value of the core radius with these 
two sets of cosmological parameters is respectively 
equal to 68kpc (ACMD model) and to 57kpc (stan- 
dard CDM model). The total binding mass de- 
rived from the X-ray data within IMpc is M tot — 
5.8 x lO 14 /ig O 1 M . This value has to be com- 
pared to those obtained from the optical gravi- 
tational lensing follow-up achieved by Fischer and 
Tyson (1997) and Sahu et al. (1998). Within 
the same radius, they provide a value of M tot — 
1.7 x 10 15 ^Mq (Fischer and Tyson 1997). 

In a previous paper, we have reported the de- 
tection of a very strong SZ effect in the direction 
of RXJ1347-1145 with the DIABOLO photometer 
(Pointecouteau et al. 1999). We have presented a 
map of the cluster central region (2' x 1'). The cor- 
responding Comptonization parameter was y(0) = 
12.7±|;? x 10~ 4 (la error bars). During the 1999 
DIABOLO run, we performed an extended map- 
ping of RXJ1347-1145 over a 4' by 4' field. The 
size of the first map was rather small so that only 
ID average slice was used to compare with the X- 
ray data. Here the size of the map allows to make 
full use of 2D information on the cluster. More- 



over, the computation of a cluster radial profile 
from a 2D map allows to discard most of the ef- 
fects induced by point sources. 

In this paper, we present the resulting extended 
map of this cluster at 2.1mm. In a first part, we 
detail the observations procedures and character- 
istics (see Sec 2). In a second time (see Sec 3), 
we describe the data processing. In section 4, we 
analyze the astrophysical data. In the last part 
(see Sec 5), we discuss the results. Throughout 
this paper, we will use the following values for 
the cosmological parameters: H = 50km/s/Mpc, 
tt m = 0.3, Q A = 0.7. 

2. Observations 

The DIABOLO instrument is a dual channel 
photometer working at 1.2 and 2.1mm. The de- 
tectors are bolometers cooled down to 0.1K us- 
ing an open cycle He- 3 He dilution refrigerator 
(Bcnoit et al. 2000). Two thermometers associ- 
ated to a heater and a PID digital control system 
are used to regulate the temperature of the 0.1K 
plate. There are three adjacent bolometers per 
channel, arranged in an equilateral triangle at the 
focus of the telescope. For a given channel, each 
bolometer is coaligned with one bolometer of the 
second channel, both looking toward the same sky 
direction. For the observations presented here, DI- 
ABOLO was installed at the focus of the IRAM 
30 meter radio telescope at Pico Veleta (Spain). 
This configuration allows us to achieve a 22" res- 
olution at 2.1mm. The 30 meter telescope focus 
being of Nasmyth type, the rotation of the field 
has to be taken into account in the sky maps re- 
construction. Desert et al. (1998) have described 
the experimental setup and reported the first DI- 
ABOLO SZ detections. 

RXJ1347-1145 has been observed in January 
1999 for a total integration time of 17 hours. 
The X-ray emission center: 0:2000 = 13 h 47 m 31 s , 
<5 20 oo = -11°45'1P (Schindler et al. 1997), has 
been taken as the map center. Each sequence of 
observation has been performed in the right as- 
cension coordinates, using the Earth rotation as 
right ascension drift such that the telescope is 
kept fixed in local coordinates. The basic obser- 
vation sequence was a 240" x 240" map in right 
ascension-declination coordinates, with a 10" dec- 
lination steps. The wobbling secondary mirror of 
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the 30m IRAM telescope has been used at a fre- 
quency of 1Hz and with a modulation amplitude of 
150" . The wobbling is horizontal (eg: at constant 
elevation), thus not aligned with the scanning di- 
rection. In order to remove systematic signal drifts 
that are produced by the antenna environment, we 
used alternatively the positive and negative beam 
to map the cluster. We have produced 89 such 
maps. 

3. Data reduction and calibration 

The data processing includes all the different 
steps described by Pointecouteau et al. (1999). 
This includes a correction from the cosmic-rays 
impacts, a synchronous demodulation of the sig- 
nal, a subtraction of the atmospheric signal corre- 
lated between the two channels, a correction from 
the atmosphere opacity, a reprojection of each sin- 
gle observation on a final RA-DEC grid, a baseline 
subtraction supported by the edges of the map and 
the computation of a final map at 1.2 and 2.1mm 
by coadding the 3 bolometers of each channel. The 
beam modulation is performed with an amplitude 
of 150" , so that some of the "OFF" positions lie 
within the limit of our map. Consequently, the 
beam switching has to be taken into account in 
the data analysis (see Sec. 4.1). 

Throughout the run, the pointing verifications 
have been performed in the direction of QSOs ly- 
ing at about the same declination as the cluster. 
The planet Mars has been used as a calibration 
target and to map the beam pattern. Assuming 
Mars is a point source with respect to the DI- 
ABOLO's beam (angular diameter around 5" in 
January 1999), the accuracy in the absolute cal- 
ibration is better than 20% at 1.2mm and 15% 
at 2.1mm. Mars has directly been mapped in 
right ascension-declination during this run. This 
method provides a direct view of the beam shape 
including eventual systematic effect as the beam 
elongation in the right-ascension direction due to 
the conjugated effects of the scanning drift speed 
and the bolometer time constant. 

4. Results and data analysis 

The reduced map for RXJ1347-1145 at 2.1mm 
is shown on Fig. la. It has been smoothed using 
a Gaussian filter with 20" FWHM. The contours 
overplotted correspond to 1, 2, 3, 4 and 5 a de- 



tection levels (negative and positive contours are 
respectively drawn with solid and dashed lines). 
The noise level on this map is lmJy/beam in a 
10" x 10" pixel (equivalent to 0.5mJy/beam in the 
22" DIABOLO beam and to 0.3mJy/bcam (equiv- 
alent to y = 3 x 10~ 5 ) in the 30" FWHM effec- 
tive beam after smoothing). This map exhibits a 
strong and extended SZ decrement in the direc- 
tion of the cluster. Unfortunately, the DIABOLO 
1.2mm data are very noisy. No positive SZ signal 
could be extract from them. We used them to sub- 
tract most of the atmospheric emission from the 
2.1mm data. 

The signal does not seem to follow a symmet- 
ric circular distribution as expected from the X- 
ray data and from the commonly used /3-model. 
The first look at the map argues in favor of sub- 
structure in the signal, and so far, in the gas dis- 
tribution. The X-ray map has been overplotted 
on the SZ map (see Fig. lc). It has been com- 
puted from the whole ROSAT/HRI observations 
(35ks exposure time) obtained from the ROSAT 
database (the map published by Schindler et al. 
(1997) included only a third of those data). 

As the SZ signal is linearly dependent of the 
gas density (scaled with n e ), the X-ray emission is 
dominated by contributions from regions of higher 
density (scaled with n 2 e ). Furthermore, the X-ray 
and SZ signals are not sensitive to the same part of 
the gas. What appears as differences between the 
SZ and the X-ray spatial distributions could be as- 
sociated to different gas phases with various phys- 
ical states (density, pressure, temperature) . In the 
following, we will compare in detail the measured 
SZ map to the expected SZ distribution deduced 
from the symmetrical cluster X-ray model. 

4.1. Modeling the SZ signal 

From the ROSAT-HRI X-ray map and the 
ASCA X-ray spectrum of RXJ1347-1145 Schindler 
et al. (1995, 1997) have extracted the clus- 
ter physical parameters: T e = 9.3 ± 0.5keV, 
n e0 = 9.4 x 10- 3 cm- 3 , 9 C = 8.4 ± 1.8" and 
/3 = 0.56±0.04 (lcr errors), assuming a /3-model for 
the gas distribution (Cavaliere and Fusco-Femiano 
1976). Using those parameters, we can compute 
the Comptonization parameter expected toward 
the cluster center: y exp {0) = (8.4±2.7) x 10~ 4 /i^ 2 . 
Following this approach, we used a /3-model to 
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describe our SZ signal. Because of the circular 
symmetry induced on the projected sky by this 
kind of model, we chose to test different /3-modcls 
on the radial profile resulting from our 2.1mm 
map (see Fig. 2). To model the SZ signal, we 
have taken into account the integration of the SZ 
spectrum on the DIABOLO passbands and the 
integration of the gas distribution over the DIA- 
BOLO beam shape. To perform a rigorous fit of 
the DIABOLO data, we need to take into account 
in the model every step of the observing proce- 
dure and of the data reduction procedure . So 
whatever the model we used, we considered it as 
a real sky and we reproduced on it all the DIA- 
BOLO's observations performed in the direction 
of RXJ1347-1147. (This includes the wobbling 
positions and amplitudes.). The resulting simu- 
lated data set of observations has been reduced in 
the same way as the real dataset by following all 
steps of the DIABOLO pipeline (demodulation, 
base line subtraction, reprojection,... see Sec 3). 

The resulting model for the flux of the SZ decre- 
ment per beam can be expressed as follow: 

F(v,H)=y(0) J TvSZ{v,T g )dv J P(Tt)L(Tt - fl')dQ' 

(1) 

where y(0) = (kBT g /m e c 2 )<7T J n e (l)dl, is the 
Comptonization parameter towards the cluster 
center (ks is the Boltzmann constant, m e the 
electron mass and c the speed of light. <tt is the 
Thomson cross section, n e (l) is the electronic den- 
sity along the line of sight). t(v) is the normalized 
DIABOLO band spectral efficiency. SZ{v,T g ) 
represents the SZ spectrum. It is a numerical 
function of v and T g with a weak dependence on 
T g that takes into account the relativistic correc- 
tions (Pointecouteau, Giard and Barret 1998). 
P( f2 ) represents the normalized gas profile pro- 
jected on the sky and L(Q) is the normalized 
DIABOLO beam shape. 

The cluster central region contains a radio point 
source, known from the NRAO VLA Sky Survey 
(NVSS) (Condon et al. 1998), which coordinates 
are a 2 ooo =13 h 47 m 30.7 s and S 2QO a =-ll°45 m 8.6 s 
(within the 3" of the X-ray center). To avoid any 
contamination and bias of the SZ signal due to 
a residual millimeter emission of this radio point 
source, we chose to include its contribution in our 



model. We considered this source as a point source 
with respect to the DIABOLO's beam. This ap- 
proach follows the method used by Pointecouteau 
et al. (1999): 

F{v,(t) = Fns{u)L(Ti) (2) 

Frs(v) is the flux of the point source. 

We have tested a total of five models on our 
data. The models A and B just include an SZ 
component. In model A, we choose the Comp- 
tonization parameter, the core radius and (3 as 
free parameters. In model B, we just let free y 
and r c and fixed (3 to 0.56 (the X-ray value). The 
three other models are a combination of an SZ 
component and a point source component. The 
parameters in model C are y, r c , (3 and Frs- In 
model D 7 we let free the y, r c and Frs param- 
eters ((3 = 0.56). Finally, for the E model, we 
have fixed the core radius and the f3 parameters 
to their respective X-ray value of 8.4" and 0.56. 
We used the y and Frs as free parameters. The 
models have been tested using a maximum like- 
lihood analysis method. The errors are obtained 
through the integration of the likelyhood function 
over the parameter space. The best fit parameters 
are gathered for the different cases in table 1. 

As shown by the results from models A and 
C, we can not constraint the core radius and the 
(3 parameters simultaneously, because those two 
parameters are too strongly coupled. For this rea- 
son (3 has been fixed to the X-ray value of 0.56 
in models B and D. Unfortunately, in those two 
cases the core radius can also not be constraint 
precisely. Its determination suffers from a degen- 
eracy and despite of the extent of our SZ map, 
we are not able to discriminate between a small 
and a high core radius model. This is due to the 
way we have defined the zero level: a baseline sub- 
traction in the RA direction fitted to 60% of the 
data point per line 30% on each edge). This oper- 
ation is needed to eliminate the low-frequency de- 
tector noises. Consequently, whatever the model 
is, the cluster extension is cut off by the baseline 
subtraction and the signal to noise ratio does not 
allow to make the difference between small and 
large core radius. For this reason, we have cho- 
sen to use the core radius value derived from the 
X-ray analysis, C — 8.4", to describe the spa- 
tial gas distribution with a /3-model. Following 
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this hypothesis, we have adopted the model E as 
the best fit model. The best fit parameters are 
?/(0) = (7.9 ± 0.5) x 10~ 4 and F RS = 5.7 ± 1.6mJy 
(with 68% confidence level error bars). In order 
to perform an accurate error analysis, we have 
also propagated the uncertainty on the X-ray pa- 
rameter determination through the whole pipeline 
of the data analysis. The extra la errors are 
Ay = 1.1 x 10~ 4 and AF RS = 1.8mJy. Because 
we just use the information on the X-ray temper- 
ature to fix the exact shape of the SZ spectra (a 
second order effect), the extra errors induced by 
the temperature uncertainty is marginally signifi- 
cant. The most important difference is driven by 
the core radius and by the [3 parameters uncer- 
tainties. Figure 2 presents the best fit radial pro- 
file model (SZ plus point source signals) overlying 
the cluster radial profile as seen by DIABOLO at 
2.1mm. The residual signal is overplotted as a 
dashed line. 

One can also perform a fit with y as a single 
free parameter. To do that, the central part of 
the map (6 C < 30") is excluded from the radial 
profile computation. The best fit parameter in 
this case is y(0) = (7.8 ± 0.4) x 10~ 4 . After- 
ward, the point source flux can be derived from 
the residual map: Fas — 5.3 ± 1.0m Jy. Those 
two values are fully consistent with the previous 
ones. (The propagation of the X-ray parameters 
uncertainties reach to respective extra errors of 
Ay = 1.0 x 10~ 4 and AF V = 1.2mJy). Our y 
value is in very good agreement within a 68% con- 
fidence level with the value expected from X-ray 
data. Moreover, it is in very good agreement with 
the value determined by Komatsu et al. (2000) 
from their 21GHz data: y(0) = (7.7 ± 1.6) x 10~ 4 . 
Finally, within the 2a of our previous determina- 
tion: 2/(0) = 12.7±|;? x 10~ 4 (Pointccouteau et al. 
1999). Combining the y(0) value expected from 
the X-ray data (y exp (0) = (8.4 ± 2.7) x 10~ 4 h 1 ^ 2 ) 
and the one we estimated from our SZ measure- 
ments, we are able to derived the H value: H = 
44 ± 6km/s/Mpc. The error quoted on H$ just 
includes the uncertainty on our y(0) determina- 
tion. If the uncertainties on the X-ray parameters 
are taken into account, the error on H become 
±15km/s/Mpc. Obviously this does not take into 
account any other uncertainty or any systematic 
errors due to the cluster geometry or to the hy- 
pothesis concerning the gas isothermality or the 



hydrostatic equilibrium. 

The residual signal resulting from the data map 
and the best fit model map difference is presented 
on figure lb (A Gaussian 20" FWHM filtering has 
been performed). The 1, 2 and 3 a detection lev- 
els have been overplotted (solid contours for the 
negative signal). This map is mostly compatible 
with noise, the higher deviation being a decrement 
(-3.1 ± 0.9mJy) located to the North-East of the 
cluster center (Aa = +35",Ac5 = +30"), with a 
significance: 3.5ct. 

Komatsu et al. (2000) have mapped RXJ1347- 
1145 at at 150GHz with the NOBA/NRO instru- 
ment. They get a 2' x 2' map of the cluster cen- 
ter. To compare our work to theirs, we have di- 
vided the center of our map in the same 4 regions 
(50" x 50" square), SE, NE, NW and SW, as they 
did. We have integrated the flux in each region. 
The value derived for each region can then be com- 
pared to the one extracted from the NOBA/NRO 
map. The results are gathered in table 2. We 
also show the integrated values for our residual 
map. We do not confirm the negative excess for 
the SE quadrant. Instead, our measurement is 
symmetrical within the noise level. If any depar- 
ture is to be searched in our map, then it is a 
positive excess in the SW quadrant: 2.8 ± 1.4. De- 
spite the low significance level of this excess and 
refering to our residual map, we can suggest that 
this excess is compatible with the presence of a 
second point source in this region. This hypoth- 
esis is supported by the point source detected at 
350GHz with SCUBA and at 8.46GHz with the 
VLA by Komatsu et al (private communication). 
Its VLA position is a 20 oo = 13' l 47 m 27.72 ;s , 6 20 oo=- 
ll°45 m 52.86 s . The SCUBA and VLA respective 
fluxes arc 17.9 ± 4.8mJy and 0.490 ± 0.043mJy. 
This source seems to be an infrared source. A mil- 
limeter residual, provided by the millimeter tail 
of a dust emission, could contribute to our sig- 
nal. Complementary millimeter observations are 
needed for this source. 

4.2. The gas mass profile 

From the DIABOLO map, we can directly esti- 
mate the gas mass which gives rise to the measured 
SZ signal. The SZ signal can directly converted 
into a gas mass value by applying the linear trans- 
formation: M(9) = C x F v {6), where C is a con- 
stant which depends on the redshift [z = 0.45) and 
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also on the gas temperature. The SZ signal map is 
obtained after the subtraction of the point source 
best fit model (see Sec. 4.1). This RXJ1347-1145 
SZ map is shown on figure lc (The X-ray con- 
tours have been overplotted) . Afterwards, it is 
converted into a gas mass map, which is then in- 
tegrated with respect to the distance to the cluster 
center (see Fig. 3). 

Due to the size of the DIABOLO map and es- 
pecially to the part used to subtract the baseline 
(30% on the edge of each line), we could only mea- 
sure the projected gas mass within a projected 
radius of 74" (ic: 0.6Mpc): M^ 8 (B < 74") = 

(1.1 ± 0.1) x lO 14 /i5 O 5/2 M . The error bar are 
quoted at a 68% confidence level. This value 
agrees very well with the estimate gas mass de- 
rived from the X-ray best fit model (under the hy- 
pothesis of a spherical j3- model): M^~ s model (r < 
0.6Mpc) = 10 14 /i- 5 / 2 M Q . 

This result confirms the agreement of the SZ 
and the X-ray data, when the hypothesis of the 
/3-model is adopted. In our case, no bias has been 
introduced in the measured projected mass by the 
baseline subtraction. In fact the part of each line 
used to subtract the baseline has been taken out 
of 0.6Mpc from the cluster center. A radius of 
0.6Mpc corresponds to an angular radius of 74.2" 
(8.8r c ). Assuming a virial radius of 10r c , the mass 
unaccounted for is estimated as no more than 16% 
of the total gas mass. 

Our determination of the gas mass agrees with 
the X-ray determination. However, the X-ray to- 
tal mass determination disagree with the strong 
lensing (Sahu et al. 1998) and the weak lens- 
ing (Fischer and Tyson 1997) determinations by 
a factor of two. Allen (1998) explained this dis- 
crepancy by the effect of the cooling flow on the 
total mass estimation through X-ray data. He con- 
cludes that the correction from this effect removes 
the differences between the estimators of the to- 
tal mass. As an example, we can propose an in- 
crease of the gas temperature by a factor of two 
(~ 18kcV). To keep the SZ signal unchanged, this 
increase has to be balanced by a decrease of the 
gas density by a factor of two. It then affects our 
results in the following ways: (a) The infered total 
mass will be divided by a factor of two (M gas (r < 
0.6Mpc) = (0.61 ± 0.04) x 10 14 h^ /2 M Q ). The 
total mass will also increase by a factor of two 



and then be in agreement with the lensing mea- 
surements, (b) As far as we suggest an increase 
of the gas temperature balanced by a decrease 
of the gas density, the SZ signal will remain un- 
changed. On the opposite, the X-ray flux will 
decrease by more than a factor of 2 (scaling the 
X-ray flux as n e (0) 2 \f[T g )). Consequently, the 
Hubble constant will be affect by a factor of 2 

(h kfI z /f x ). 

The previous scheme is very simple. The 
thermodynamic state of the intracluster medium 
is probably more complex mainly made of two 
gas components. A cold (oc 9kcV), dense and 
peacked central component surrounded by a hot 
(oc 18kcV), thin and extended component. In this 
case the cluster should present a non isothermal 
radial profile (assuming the spherical symmetry). 
This kind of hypothesis could explain the very 
strong degeneracy which we encountered while 
trying to determine the core radius (see Sec. 4.1). 
Unfortunately, the quality of our data is not good 
enough to argue in this direction. Nevertheless, 
this suggestion is not unlikely. Recent XMM- 
Newton observations have already shown such a 
behavior in other cooling flow clusters: A1795 
(Arnaud et al. 2000; Tamura et al. 2001), A1835 
(Peterson et al. 2000) 

4.3. Point sources 

The best fit model gives a flux of 5.7 ± 
1.5mJy/beam for the central point source. This 
point source can be visualized on the data map 
after the subtraction of the SZ component. Fig- 
ure Id shows clearly an excess of positive emis- 
sion in its center. The NVSS map (Condon et 
al. 1998) has been overplotted as contours. Both 
signals match very well in term of position. 

This point source has already been observed at 
different frequencies. It has been detected in the 
NVSS by Condon et al. (1998) with a flux of 
46.9 ± 1.9mJy at 1,4GHz. Komatsu et al. (1999) 
have reported various measurements of this point 
source: an OVRO flux of 10.1 ± 0.5 at 28.5GHZ, 
a NMA (Nobeyama Millimeter Array) flux of 5 ± 
1.5mJy at 100GHz and a 2a upper limit of 4.8mJy 
obtained with SCUBA at 350GHz. From this set 
of measurements and assuming a power law spec- 
trum for the point source radio emission, they 
have derived: F RS = (55.7 ± 1.0)^f ±0 02 mJy. 
More recently, Komatsu et al. (2000) have ob- 
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served it with the VLA at 8.46 and 22.46GHz. 
The respective detected flux are 22.42 ± 0.04 and 
11.5 ± 0.17mJy. 

Using those latest measurements, we update 
the point source power spectrum determination: 
Frs = (77.8 ± l.7)vQ° H f ±0 ' 01 mJy. The qual- 
ity of this fit is moderate, due to the very ro- 
bust constraints given by some of the data points 
(see Fig. 4). It is difficult to conclude firmly 
that a single power spectrum can model correctly 
the synchrotron emission from radio wavelengths 
down to millimeter wavelengths. The intensity 
and the spectral shape of the synchrotron emission 
is driven by the internal magnetic field of the radio 
source. Furthermore its strength also drives the 
maximum frequency of emission allowed. To date, 
a lot of radio sources have been observed with var- 
ious intensity and spectral shapes. Because of this 
high rate of variation from one source to an other 
and because of the lack of millimeter data, we are 
not able to state on an eventual cut-off in the spec- 
tral shape of our radio source. For this reason, we 
adopted the previously defined power law model. 
It provides us an estimation of the flux contributed 
by the central point source at 143GHz (2.1mm). 

To explain the nature of the millimeter emis- 
sion we have detected, we will now compare the 
residual radio flux expected at 143GHz (2.1mm) 
and 350GHz (850/xm) from the previous power law 
spectrum to the flux respectively obtained from 
the DIABOLO data and from the SCUBA data. 
At 2.1mm, we have extrapolated: (2.1mm) = 
4.3±0.3mJy. This estimation is at la compatible 
with our determination from the DIABOLO map. 
The 350GHz (eg: 805^tm) flux can be estimated 
too: (0.85mm) = 2.7 ± 0.4mJy. Now, from 
the SCUBA data points published by Komatsu et 
al. (1999), we have subtracted the positive SZ 
contribution, computed from the cluster SZ spec- 
trum and scaled by our y best fit value. The SZ 
radial profile has been obtained from the model C 
(see Sec. 4.1) convolved with the SCUBA beam (a 
Gaussian beam with a fwhm = 15"). We fixed 
the value of the signal offset (so-called DC off- 
set) to 2. 7m Jy /beam as published by Komatsu 
et al. We finally deduced for the point source: 
F RS (0.85mm) = 1.8 ± l.OmJy. The low signal to 
noise ratio of this last flux does not allow us to 
consider this result as the detection of a real sub- 
millimeter point source. 



Nevertheless, our DIABOLO flux and the 
SCUBA flux are both compatible with the esti- 
mation extrapolated from the point source power 
spectrum. For this reason, it seems reasonable to 
deduce that the point source emission seen in our 
map is likely due to the synchrotron millimeter 
tail of the central radio point source emission. 

5. Conclusion 

We have produced an extended map (4" x 4") 
of the RXJ1347-1145 cluster with the DIABOLO 
photometer. We have drawn the distribution of 
the SZ signal up to 74" (0.6Mpc) from the cluster 
center. This SZ signal is as much extended as 
the X-ray emission and is the strongest SZ effect 
detected to date. 

The SZ map allows us to directly derived the 
projected mass distribution and by the way to 
measure the cluster gas mass up to 0.6Mpc. We 
derived M gas = (1.1 ± 0.1) x 10 14 h^ /2 M Q . The 
gas mass estimated within the same region, un- 
der the assumption of a classical /3-model and the 
cluster parameters derived from the X-ray data, 
agrees with our value. 

Due to the map noise level and to the observing 
strategy, we are not able to firmly conclude to the 
existence of substructure in the gas distribution. 
Nevertheless, our map presents some asymmetric 
features at more than a 2a level. 

At the moment, the first results of the Chandra 
and the XMM-Newton satellites, concerning the 
observation of galaxy clusters, already show some 
departure in the density distribution from the j3- 
model symmetry, as well as in the temperature 
distribution from the isothermality (Fabian et al. 
2000; Vikhlinin, Markevitch and Murray 2000; 
Arnaud et al. 2000; Tamura et al. 2001; Peterson 
et al. 2000). For this reason, the comparison of 
the upcoming Chandra and XMM-Newton X-ray 
obervations to the actual and the future SZ data 
(from interferometers and bolometer arrays) are 
needed to understand the cluster physical struc- 
ture. 
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Fig. 1.— (a) 2.1mm map of RXJ1347-1145 ob- 
tained with the DIABOLO photometer. This map 
and the three others have been smoothed by a 
gaussian filter with 20" FWHM. The la noise level 
in term of flux is lmjy in a pixel of 10" x 10" 
(equivalent to 0.5mJy/beam in the 22' DIABOLO 
beam and to 0.3mJy/beam in the 30" FWHM ef- 
fective beam after smoothing), (b) Residual map 
after subtraction of the best fit model (SZ plus 
point source, see text). In the maps a and b the 
overlying contours correspond to the 1, 2, 3, 4 
and 5 a detection levels. The negative and pos- 
itive contours are respectively drawn with solid 
and dashed lines, (c) RXJ1347-1145 SZ signal (af- 
ter subtraction of the point source best fit model) 
compared to the X-ray cluster emission. The dif- 
ferent contours overplotted correspond to 3, 5, 10, 
30, 50, 70 and 90% of the X-ray maximum of emis- 
sion, (d) Point source detection at 2.1mm after 
subtraction of the SZ best fit model. The con- 
tours overlying correspond to the NVSS map (1, 
10, 30, 50, 70 and 90% of the radio maximum 
emission). The central position in the four map 
correspond to the position of the X-ray center: 
a 20 oo = 13 h 47 m 31 s , 5 20Q0 = -ll°45'li". 

Fig. 2.— RXJ1347-1145 radial profile at 2.1mm 
computed from the DIABOLO 2.1mm map. The 
data point are plotted with their associated la er- 
ror bars. The best fit model (including an SZ com- 
ponent and a point source component, see text) is 
overplotted (solid line) . The residual radial profile 
is plotted as a dashed line. 

Fig. 3. — RXJ1347-1145 gas mass profile com- 
puted directly from the 2.1mm map. The solid 
line represents the gas mass profile of the SZ best 
fit model. The integral of the DIABOLO beam 
(dashed line) has been overplotted and scalled to 
the gas mass profile to show the extension of the 
SZ signal. 

Fig. 4. — Central point source spectum from 1.4 
to 350GHz. 
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Table 1: Best fit parameters. 



Models y(0) [1CT 4 ] 6 C (arcscc) f3 F RS (mJy) xV(" - l T 



A 


4 1 +i - Y 


56.3™ 


u - oy -0.61 


b 


1.1 


B 


c Q+0.8 
°-°-0.3 


32-41^ 


0.56 b 


b 


0.9 


C 


5-7±?;g 




n 5i +0 - 34 

u - 01 -o.n 




1.1 


D 




14.8±f 7 7 


0.56 b 


LI l'^ 


0.9 


E 


7-91S1 


8.4 b 


0.56 b 




1.1 



"where (n-1) is the number of degrees of freedom 
6 fixcd parameter, see text 
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Table 2: Fluxes in 50" x 50" SE, NE, NW and SW 
regions as denned by Komatsu et al. (2000). The 
column respectively gather the DIABOLO fluxes 
(1), the residual flux after the subtraction of the 
best fit model (2), the NOBA/NRO fluxes (3). 
Columns (1) and (3) have to be compared. The 
second column try to highlight some eventual ex- 
cess (positive or negative). For each collumn the 
associated la noise flux is quoted. 



Regions 


map 


residual map 


NOBA/NRO 




(mJy) 


(mJy) 


(mJy) 


SE 


-12.2 


-0.8 


-11.3 


NE 


-9.9 


-0.7 


-4.7 


NW 


-10.1 


-0.4 


-3.3 


SW 


-6.8 


2.8 


-6.1 


la noise 


1.4 


1.4 


2.0 
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13 




R (Mpc) 



14 



100 



+ 



10 



* SCUBA 

A DIABOLO 

□ NMA 

O OVRO 

X VLA 22.46Ghz 

O VLA 8.46Ghz 

+ VLA 1.4Ghz 
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100 
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